APS/123-QED 



Enhancement of Nuclear Spin-Lattice Relaxation Rate and Spin Susceptibility 

due to Valence Fluctuations 
-Origin of Anomalously Enhanced Wilson Ratio in Ce and Yb Systems- 



o 
o 

(N 
(N 



I 

o 
o 



00 

m 
o 

On 
O 



Shinji Watanabe^ and Kazumasa Miyake^ 
Department of Applied Physics, University of Tokyo, 
Hongo 7-3-1, Bunkyo-ku, Tokyo, 113-8656, Japan} 
Division of Materials Physics, Department of Materials Engineering Science, 
Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japar? 

(Dated: November 24, 2008) 

We show theoretically that the nuclear spin-lattice relaxation rate (TiT)^^ and the spin suscep- 
tibility XsiT) exhibit divergent behaviors toward zero temperature at the quantum critical point 
(QCP) of the first-order valence transition. Remarkable enhancement in (TiT)"^ and Xs{T) is in- 
duced by valence fluctuations even at the valence-crossover temperature far away from the QCP. 
This mechanism well explains peculiar behaviors observed recently in YbAuCu4 and also gives a sys- 
tematic explanation for YbXCu4 for X=In, Au, Ag, Tl, and Pd from the viewpoint of the closeness 
to the QCP. This also explains anomalously enhanced Wilson ratio observed in the paramagnetic Ce 
and Yb based compounds. This offers a new concept that spin fluctuations are induced via relative 
charge fluctuations, which can be generally applied to the systems with valence instabilities. 
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Quantum critical phenomena in strongly-correlated 
electron systems have been discussed extensively in the 
context of magnetic phase transitions [ij, ^2, ^] . Recently, 
instabilities in charge sectors have attracted much atten- 
tion, since underlying influence of valence instability is 
suggested by variety of materials [4, ,5;, §, 7] : Importance 
of critical valence fluctuations has been argued as a pos- 
sible origin of anomalies such as T-linear resistivity, en- 
hanced residual resistivity and superconductivity for the 
materials vi^ith valence-fluctuating ions such as Ce and 
Yb. 

Valence transition is isostructural transition as known 
as 7-q; transition in Ce metal [8] and also in YbInCu4 [9j 
characterized by a jump of the valence of the Ce and 
Yb ion. In YbInCu4 the first-order valence transition at 
T ~ A2K takes place with the valence of Yb being -t-2.97 
(-1-2.84) in the high (low)-temperature phase Inter- 
estingly, anomalous behaviors in the YbXCu4 have been 
revealed by systematic measurements: A remarkable 
enhancement in the nuclear spin-lattice relaxation rate 
(TiT)-! has been discovered recently for X=Au 
A mysterious behavior of this material has been recog- 
nized in the spin susceptibility, which increases enor- 
mously toward zero temperature [isj in spite of the 
antiferromagnetically-ordered ground state with Tn ~ 
0.8 K By applying the magnetic field Tn is sup- 

pressed to K at iJc 1-3 T and a remarkable point is 
that by further applying H the enhancement of (TiT)~^ 
emerges at finite temperature, T — T^{H) Since 
these anomalies appear even far away from He, i.e., in the 
regime where the magnetic order is completely destroyed, 
this is not due to the magnetic fluctuations. 

Furthermore, it has been discovered that the ^^Cu 
NQR frequency z^q shows a sharp drop at when T 
decreases ^Ij- Since i^q measures the charge distribu- 
tion of the Yb and surrounding ions, the change of vq 



indicates that the Yb valence changes at sharply. A 
peak structure in the spin susceptibility Xa{T) at Tv has 
been also detected for X=Ag [13i] and Tl [iJl, and en- 
hancement in (TiT)~^ has been also found at Tv for 
X=Ag jT^ and X=Pd recently. These observations 
suggest that this behavior is not specific to a special ma- 
terial, but is rather universal. Since valence fluctuations 
are ascribed to the relative charge fluctuations between 
f and conduction electrons, these observations challenge 
the conventional concept that the nuclear spin-lattice re- 
laxation rate and the spin susceptibility reflect magnetic 
fluctuations in the paramagnetic-metal phase. 

In this Letter, we resolve this puzzle by showing that 
valence fluctuations indeed induce spin fluctuations. We 
show that (TiT)^^ as well as Xs{T) shows divergence 
toward zero temperature at the quantum critical point 
(QCP) of the valence transition. Even in the valence- 
crossover region at finite temperatures away from the 
QCP, (TiT)~^ and Xs(T) are shown to be enhanced. 
This mechanism gives a systematic explanation for pe- 
culiar behaviors observed in YbXCu4 for X=In, Au, Ag, 
Tl and Pd, and also accounts for anomalously enhanced 
Wilson ratio observed in several Yb and Ce materials lo- 
cated near the QCP of the valence transition. This result 
offers a new concept that spin fluctuations are induced 
via relative charge fluctuations, which can be generally 
applied to the systems with valence instabilities. 

Let us start our analysis by introducing a minimal 
model which describes the essential part of the Yb and 
Ce systems in the standard notation [l^, [l^ : 



H — He + Hf + Hhyh + Hu{^ 1 

Ek^Ek^aCka, Hi ^ EiY. 
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Us J2f=i nl^^li^ ^hyb = Eka (/Lck^ + cl^fka) and 

The Ufc term is the Coulomb 
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repulsion between f and conduction electrons, which is 
considered to play an important role in the valence tran- 
sition: In the case of Ce metal which exhibits the 7- 
a transition, the 4f- and 5d-electron bands are located 
at the Fermi level [l^l- Since both the orbitals are lo- 
cated on the same Ce site, this term cannot be neglected. 
For YbInCu4, the considerable magnitude of the In 5p 
and Yb 4f hybridization was pointed out by the band- 
structure calculation [2l| and recent high-resolution pho- 
toemission spectra has detected a remarkable increase 
of the p-f hybridization at the first-order valence tran- 
sition [22|. These results suggest importance of Vk and 
Ufc- Actually, the reason why the critical-end tempera- 
ture is so high as much as 600 K in Ce metal in contrast 
to that in YbInCu4 can be understood in terms of Ufc'. 
In YbInCu4 , Ufc is the intersite interaction, which should 
be smaller than that of Ce metal. This view also gives 
an explanation for the reason why most of Ce and Yb 
compounds only shows the valence crossover. Namely, 
most of the compounds seems to have a moderate value 
of Ufc due to its intersite origin, which is smaller than the 
critical value to cause a jump of the valence. However, 
even in the valence-crossover regime, underlying influence 
of the valence instability causes intriguing phenomena as 
shown below. It is noted that importance of Ufc has been 
discussed by several authors for YbInCu4 [H, [IJ] . 

In this model ^ the first-order transition between the 
larger (uf) and the smaller (nt) is caused by Ufc, since 
the large Ufc forces electrons to pour into either the f 
level or the conduction band 0, [13] ■ Figure [TJa) shows 
the ground-state phase diagram in the Ef-t/fc plane de- 
termined by the density-matrix renormalization group 
(DMRG) for ek = -2cos(A:), V = = 0.1 and U ^ 100 
at filling E»=i("« +0/(2A^) = 7/8 jj. We note that 
essentially the same phase diagram has been obtained 
in the infinite dimensional system [26| . The first-order- 
transition line (brown line) separates the paramagnetic- 
metal phase, namely, the larger (rif ) phase in the smaller 
Ef and Ufc regime, and the smaller (uf ) phase in the larger 
Ef and Ufc regime. As Ufc decreases the jump in Uf at 
the first-order transition decreases and terminates at the 
QCP, at which the valence fluctuation diverges. 

The valence fluctuation is measured by the dynam- 
ical valence susceptibility defined by x^i^li'i^n) = 
/(f dr(T.,nf(q,r)nf(-q,0))e*""'" with u;„ = 2mTT {n = 
0, ±1, ±2, • • • ). The most dominant part of near the 
QCP is expressed by the four-point vertex function F 
in Fig. mja), which satisfies the integral equation shown 
in Fig. [DJb). Here, the solid line and dashed line rep- 
resent quasiparticle parts of Green functions near the 
Fermi level G'^(q,ie„) ~ l/(ie„ - S*^), Gf{q,ien) 

— ^-qcr); respectively, where £„ = (2n -I- 1)ttT 
and the origin of the energy is set as n = 0. Here, 
satisfies {z - Sf - - Eq) - | = 

{z - where ^(z) = + with S^, 

and Sqg. being the self energies due to Uff and Ufc, 
respectively. The effect of the renormalization ampli- 
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FIG. 1: (Color) (a) Ground-state phase diagram determined 
by the DMRG for = -2 cos(it), V = = 0.1 and U = 100 
at filling E»^i(Wi + ni)/{2N) = 7/8 7]. (b) Valence suscep- 
tibility has a peak on the dashed line in (a), (c) System-size 
dependence of the total-charge compressibility at the QCP 
{Ei,Ufc) = (—4.5206,5.8460). (d) Schematic phase diagram 
in the Sf-Utc-T space for a certain Vk and Uh{> Utc). The 
first-order valence transition surface (brown surface) with the 
critical end line (blue line) touched on T = at the QCP con- 
tinues to the valence-crossover surface (light-green surface). 
The valence-crossover surface at which k\T) has a maximum 
is denoted as Tv (see text). 

tude defined by a^^ = [1 - dJ:lJe)/ds + |Fq^(e)|V(e - 
e^)'-2Re[V*A£)d^U^)/de]/{e-e^)r^U^o and a- = 
Cq^lVq^ (0)p/eq is absorbcd in the definition of vertex 
functions, internal and external [IJ. The vertex F con- 
sists of the interactions via Ufc with vertex parts includ- 
ing Uff, whose contributions are represented by the dou- 
ble wigly line as Ufc in Fig. [2^c) and all the other con- 
tributions such as entangled Ufc and Uff shown as Fjl^ 

in Fig. Hfd). Here, f = a[/fc/x", where a and are 
the 2x2 matrices with elements [a]a-a' = a^'^Sv and 
[x"](T(t' = Xcrcr' defined as Fig[2Kd), respectively. / is the 
unit matrix. Then, we have oc Tr[%"(/ — C/fc/x")""'^]- 
By using the relation Xchargc = XtT + XU and xi'pin = 
XTT - Xn and noting the fact that Xchargo < xi'pin 
holds for typical heavy-electron systems |2^], we have 

X*^of xipi„/(l-^fcxipi„)■ 
Figure md) illustrates the schematic phase diagram 
in the T-Sf-Ufc space. At the critical end line of 
the first-order transition illustrated by the blue line in 
Fig. [IJd) where the denominator of equals zero, 

1 — ^fcxipin(0' 0) = 0' the valence susceptibifity = 
X^{0, 0) = —d{n^)/d£f diverges. Note here that shows 
not only divergence at the critical-end line, but also has 
a maximum at the valence-crossover surface (k^ has a 
peak in Fig. [IJb) on the dashed line in Fig. [H^a)), which 
is illustrated as the surface extended from the first-order 
one in Fig. [Tfd) 0, HI]. This implies that even if the 
system is away from the critical end line, the valence 
fluctuation is enhanced at the crossover surface Tv in the 
parameter space of T and pressure and/or chemical dop- 
ing. Furthermore, it has been revealed recently that even 
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FIG. 2: (a) The most dominant part near the QCP of and 
X+- with four-point vertex F and (b) its integral equation, 
(c) Vertex parts including Ug (hatched area) connected by Utc 
(wiggly line). All contribution of this type is represented by 
the double wiggly line as C/fc. (d) F for a^nd (e) for X+-- 
is antiparallel spin to a. Fh contains all contributions other 
than the diagrams expressed in (c) . (f ) RPA-type self energy 
for Ucc- The solid line and dashed line represent and G^^^, 
respectively. 



in the valence-crossover regime, the critical point is in- 
duced by applying the magnetic field poj . This explains 
the role of the magnetic field in the T-H phase diagram 
of YbAuCu4 ^ and YbPdCu4 [13 ■ Namely, the mag- 
netic field not only destroys the magnetic order, but also 
plays a role to make the valence-crossover temperature 
Tv finite. This is consistent with the experimental fact 
that Ty{H) emerges under H even much larger than He 
in the T-H plane, at which vq shows a sharp change. 

Now let us focus on the nuclear spin-lattice relax- 
ation rate defined by ^ = Ea |i„p 

where 7n is the gyromagnetic ratio of nuclear spin, 
gi is a Lande's g factor for f electrons and Aq 
is the hyperfine-coupling constant. The dynami- 
cal f-spin susceptibility is defined by x+-{^j^'^n) = 

of _^ of / „ nWai^^n 



dr(T^S'^(q,r)S'L(-q,0))e*""^. We note here that 
X+- fia-s essentially the same structure as x^: whose 
most dominant terms near the QCP is expressed in 
Figs. [2Ja) and (b). The vertex f in Fig. ^h) is given 
by f = a^ff t^fcXipin as shown in Figs. M.e) and (c). 
Then, we have xl^ « XspJi^ ~ UtcXspin)- Since SU(2) 
symmetry of the system ensures xipin — xi'pinj turns 
out that x+- fia-s the same form as ■ Hence, at the 
critical-end line as well as the QCP of the valence tran- 
sition (blue line in Fig [Hd)), 1 - f/fcxipi„(0, 0) = 

holds, which makes x+-(O:0) diverge |30j. Since the 
spectral weight of f electrons is dominated by the inco- 
herent part around e ~ et, the coherent part responsi- 
ble for quasiparticles amounts to a tiny contribution in 
the order of a^^^ ~ "iband/^i* ^ 1- Then, the q de- 
pendence of x^ijj(q, zcjo) is estimated as Xj;|pijj(q, iwo) ~ 
Xspin{0,0) [l + Aq'^ -iCujo/q] with A in the order of 
0{Vkp /leil)"^ < 10~^ for typical heavy-electron systems 



with the spherical Fermi surface with q = |q| [J]. 
Then, the most singular term is evaluated to show 
the (TiT)^^ ~ q'^~^ divergence for g ^ in the d- 
dimensional system. This result is in sharp contrast to 
the single-orbital system where charge instability does 
not occur simultaneously with the spin instability. 

The uniform spin susceptibility is defined by Xs{T) = 
xi{T) + XsiT) with Xs = dm^/dH\H=o and = 
Y^i{Sf^)/N, when the magnetic field is applied to 
(P) as -5fMBi?E,5p - gcfiBHJ2^S'i^. By us- 
ing the dynamical spin susceptibility x+'-ili'^^n) = 

/(f rfT(T^S'^(q,T)S'^(-q,0))e*'^"^, Xs{T) is expressed as 

Xs{T) = |m| [ghl- + 9i9.X+- + 9c9txt- + 9lx+-] , 
where q and w„ are set to be zero in the right hand 
side (r.h.s.). In heavy-electron systems, the uniform 
spin susceptibility is dominated by the f-electron part 
as Xs(0) » Xs(0) M- Then, gfx%- + 9i9cX+- » 
9c9fX+- + ffcX+- holds. Since x^- = xf-, we have 
Xs{T) - {3/2)nlgfxl-- Hence, at the QCP, Xs(0) shows 
uj~^ divergence with ujy = 1 - i/fcxipi„(0, 0). 

We note here that quite different behavior appears in 
the charge sector: Our DMRG calculation apphed to the 
Hamiltonian ([T]) shows that even at the QCP of the va- 
lence transition where diverges, the total charge com- 
pressibility Ktot = d{n^) /dfi + d{n'^)/diJ, does not diverge 
(see Fig.[Tl[c)). Namely, the r.h.s. of /«tot(0) = X^ + X^'^ + 
-^cf _|_ ^cc g^j.g cancelled each other to give a finite value in 
spite that x^ diverges. Here, q and cun are set to be zero 
in x"''(q,ic^„) = /„''dT(r,n^(q,T)n''(-q,0))e'""^ This 
is in sharp contrast to the mean- field result, where the 
first-order valence transition is accompanied by the phase 
separation. Namely, diverging relative charge fluctuation 
(i.e., valence fluctuation) also induces the instability of 
the total charge in the mean-field framework. However, 
our finding shows that quantum fluctuations and elec- 
tron correlations can make the total charge stable even 
though the relative charge unstable. This is ascribed to 
the fact that the order parameter of the valence tran- 
sition, — Y^'iLi is i^ot the conserving quantity, i.e., 
[n^ ,H]^ ^ [7]. This result indicates that the first-order 
valence transition is not accompanied by the phase sepa- 
ration at least in electronic origin. We point out that the 
g factors in XsiT), which usually differs between f and 
conduction electrons, i.e., gi ^ gc prevent from the can- 
cellation among x+'- as occurred in the r.h.s. of Ktot(O), 
and hence the divergent behavior can emerge in Xs(0) in 
contrast to Ktot(O). 

The enhanced Xs {T) toward low temperature has been 
observed in YbAuCu4 The NQR and susceptibility 
measurements suggest that YbAuCu4 is located near the 
QCP of the valence transition at H = although 
the antiferromagnetic order masks the ground state with 
Tn = 0.8 K. Namely, the valence-crossover temperature 
Tv is suppressed to be close to T = K so that Xs{T) 
is interpreted to be enhanced toward low temperature 
by the critical valence fluctuations. On the other hand, 
in YbAgCu4 where no magnetic transition has been ob- 
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served at ambient pressure at = 0, the lattice con- 
stant changes around = 40 K, suggesting the valence 
crossover |15|. A remarkable point is that Xs{T) has a 
peak at Tv |13l |. whose maximum value is one order of 
magnitude smaller than that of YbAuCu4. This is again 
consistent with our theory, since at finite Tv, i.e., at the 
valence-crossover point Xs(T'v) is enhanced but does not 
diverge, while in the case of Tv = 0, i.e., at the QCP 
Xs(0) diverges (see Fig. [TJb) and Fig. [Hd)). We note 
that a recent {TiT)~^ measurement has also detected a 
broad peak around in YbAgCu4 [3 ■ It is noted that 
the maximum in XsiT) has been also observed at Tv for 
X=T1 [l3|. Hence, as X moves as An, Ag to Tl, Tv in- 
creases, i.e., the distance from the QCP becomes long, 
which makes the peak value of the spin susceptibility 
Xs(Tv) small. Then, our theory explains these system- 
atic observations quite consistently. 

Finally, we argue the Wilson ratio near the 
QCP. The self energy for f electrons is given 

hv (n ip — H 2l7fcG™(p-q,te„-iaJn.) 

^o-lP, - JV Z^q.m l-(7toX?<,(q,i'^™) 

within the RPA as illustrated in Fig. [Iff). Here, 
Xaai^^i^m) = - ^ Ek,n (k, ie„)G^'= (k q, ie„ iw™) 
where G^^{k,ien) = - + and 

G-(k,ze„) = - + with C/fc 

set to be zero in the previous definition of -Bqg.- 
Here, ^ [1 - di:l^{e) / de\,=o + ^kV4]"' and 

<<T = a^^ykV^L and G?^ and Gf^^ denote the 
incoherent parts. Near the QCP small q and uj com- 
ponents are important and the denominator of is 
expanded as 1 — UicX^^a[(\^ ^ + ^^Xt" "^1^ ~ iCuu/q 
with cjv = 1 - J7fcxi'=„(O,0) g. Then, we eval- 
uate r7i*/"iband = 1 - 9ReS^(pF,£)/i9e|£=o = 
1 - (a^;,i7fc/(4^2^F^))ln|(iv/(t^v + Aq,'')\ to the 
leading order of a!^^ ~ ^pfo-' where v-p is the Fermi 
velocity and qc is a cut-off. Namely, the Sommer- 
feld constant 7c shows a log divergence at the QCP. 
Since Xs(0) shows a uj~^ divergence, the Wilson ratio 
i?-vv — 47r^fc|xs/(3(5fAfB)^7c) diverges at the QCP of the 
valence transition (<Dv — s- 0). 



Although the Gaussian fixed point of the QCP of the 
valence transition ensures the validity of the RPA de- 
scription [1] , more qualitatively (Dv should be determined 
self-consistently as done in the SCR theory for spin fluc- 
tuations pj. This kind of analysis can be executed for 
valence fluctuations starting from the valence suscepti- 
bility x*(q, w)""'^ (Dv -I- Aq^ — iCqUi, expanded near 
(q, = (0,0). Here, Gg is given, in the T ^ hmit, by 
the form C/max{q,l~^} with I being the mean free path 
of the impurity scattering [s^l . Then, in the realistic sit- 
uation, the dynamical exponent Zd = 3 is expected to be 
observed at low T, except in the very vicinity of T = K 
where impurity scattering is dominant and hence the sys- 
tem is described by Zd = 2 [s^. Since in the three di- 
mension with Zd = 3 the susceptibility and Sommerfeld 
constant behave as Xs{T) ^-4/3 g^j-^^j ^ — InT, re- 
spectively, in case that incoherent part of f electrons gives 
minor contributions to the criticality, the Wilson ratio di- 
verges in this case, too. Namely, our results based on the 
RPA are considered to be qualitatively correct and more 
quantitative arguments for comparison with experiments 
will be reported in the separated paper [s^j- 

We note that Xs(0) and 7o extrapolated to T ^ K 
from the T > Tn = 0.8 K data in YbAuCu4 shows i?w ^ 
4.5 [l^, which exceeds Rw — 2. Enhanced Wilson ratio 
i?w ^ 3 has been also observed in Ceo.9-a;Laa;Tho.i with 
a; = 0.1 [m where the 7-Q;-transition temperature is sup- 
pressed closely to zero temperature. These materials are 
located near the QCP of the valence transition and hence 
our theory gives an explanation for these enhancements. 
Enhanced i?w has been also observed in other para- 
magnetic materials such as YbRh2(Sio. 95060.05)2 [ssl ]. 
YbIr2Si2 [M] and /3-YbAlB4 [37], which also suggests un- 
derlying influence of valence fluctuations. 

We stress that our results can be generally applied to 
the systems with valence instabilities. Experimental ex- 
amination of our predictions is highly desired. 

The authors thank S. Wada and S. Nakatsuji for show- 
ing us their experimental data prior to publication. 
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